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1^ ' 

The low-energy electronic structure of icosahedral fullerenes is studied within the field-theory 
, model. In the field model, the pentagonal rings in the fuUerene are simulated by two kinds of gauge 

fields. The first one, non-abelian field, follows from so-called K spin rotation invariance for the 
spinor field while the second one describes the elastic flow due to pentagonal apical disclinations. 
For fuUerene molecule, these fluxes are taken into account by introducing an effective fleld due to 
magnetic monopole placed at the center of a sphere. Additionally, the spherical geometry of the 
04 , fuUerene is incorporated via the spin connection term. The exact analytical solution of the problem 

(both for the eigenf unctions and the energy spectrum) is found. 
PACS: 73.22.-f 
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I. INTRODUCTION 



The electronic structure and elementary excitations of 
fuUerene molecules have been of steady interest since the discovery of the first so-called buckminsterfuUcrenc Cgo [1] 
because the knowledge of the electronic structure gives an important information about the electric and photo con- 
C 1 ductivity, magnetic behavior, etc. All these characteristics were found to be rather unique in fuUerene molecules (see, 
' e.g., review [2]), which can be effectively used in some practical applications in devices based on fullerenes. 
: There is a number of different theories to study this problem, which can be roughly divided into three general 
O groups. The first one includes empirical methods like the free-electron gas [3], and tight-binding [4~7] calculations. 
^ O^ . The second one uses ab initio quantum chemistry calculations [8]. The third group considers the continuum models 

within the effective- mass description [9-11]. 
ps) I While the continuum description is limited to the electronic states close to the Fermi level, it has some interesting 
^ ■ attractive features. First of all it gives a possibility to study large fullerenes where the numerical analysis is a rather 
^^O difficult task. Second, it reveals the long-distance physics which is of importance in various carbon nanoparticles. 
Finally, the continuum description allows to elucidate "true" topological effects like the appearance of the Aharonov- 
Bohm phase and anomalous Landau levels due to disclinations. 

In this paper, we formulate a continuum model to study low-energy electron states in icosahedral fullerenes. The 
ly-j ' model is a variant of the effective field theory on a sphere describing Dirac wavefunctions interacting with two types of 
gauge fiuxes. One of the fluxes is due to so-called K spin rotation invariance (see [12] for details) and the second one 
"Y^ comes from the local S0(2) invariance of the two-dimensional elastic Lagrangian in the presence of disclinations [13]. 
Ci Actually, the second flux describes the elastic flow through a surface due to a disclination and has a topological origin 
^ (its circulation is determined by the Frank index, the topological characteristic of the defect). For this reason, this 
I flux exists even within the so-called " inextensional" limit (which is usually adjusted to fuUerene molecule [14]). Notice 
that the topological origin the elastic flux results in appearance of the disclination-induced Aharonov-Bohm-like phase 
(see [15]). 

^ , It should be mentioned that the first continuum model for the description of fuUerene molecule was presented 
in two papers [9,10]. A different variant of the continuum model for fuUerene was suggested in [11]. The first 
model neglected the topological origin of the disclination defects while the second model missed the K spin rotation 
I invariance. However, in both papers the non-trivial zero- mode electronic states are considered. In this paper we found 
5— I ' an exact analytical solution of the problem at low energies. 

■ 



II. GENERAL FORMALISM 



Let us start from the standard formalism based on the effective-mass theory proposed in [16] to study the screening 
of a single intercalant within a graphite host. A graphite host is considered as a single graphite plane (graphene). 
Actually, the effective-mass expansion is equivalent to the k ■ p expansion of the graphite energy bands about the K 
point in the Brillouin zone when the intercalant potential is zero. In fact, there are two degenerate Bloch eigenstates 
at AT, so that the microscopic wave function can be approximated by 
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(1) 



where k = K + k'. Keeping terms of order k' in the Schrodinger equation one can obtain the secular equation for 
functions f\fi{k') and after diagonaUzation one finally gets the two-dimensional Dirac equation (see [16] for details) 

-iai'd^iPif) = Ei,{r) (2) 

where tr'* are the conventional Pauli matrices (/i = 1, 2), the energy E is accounted from the Fermi energy, the Fermi 
velocity Vp is taken to be one, and the two-component wave function xp represents two graphite sublattices {A and 
B). As was mentioned in [16] the essence of the k- p approximation is to replace the graphite bands by conical 
dispersions at the Fermi energy. In addition, one has to take into account two independent wave vectors {K^ = K 
and K- = —K) in the carbon lattice which give the same conical dispersion. Therefore, the states | K±A > and 
I K±B > can be chosen as the full basis set [12]. 

For our purpose, to take into account both the spherical geometry of a fuUerene molecule and disclinations on its 
surface we have to modify the model (2) by introducing three additional fields. 

A. The compensating fields 

As a first step, let us consider a single disclination in a graphene sheet. To introduce the fivefold in the hexagonal 
lattice one has to cut the sector 27r/6 and glue the opposite sides {OB and OB' in Fig.l). 




FIG. 1. A disclination in elastic membrane (right) and the image of the membrane (left). 
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Notice that this is a typical cut-and-glue procedure to incorporate a discUnation in an elastic plane. There exist several 
continmim theories for the description of disclinations in flexible membranes (sec [f3,17,f8]) which allow to formulate 
similar von Karman equations. In the " inextensional" limit (free to bend, impossible to stretch), the exact solution for 
an isolated positive disclination was found in [17]. It was shown that the elastic membrane becomes buckled and takes 
a conical form (see Fig.l). For our purpose, the most appropriate is the gauge theory of disclinations on fluctuating 
elastic surfaces formulated in [13]. Indeed, this is the most convenient theory to take into account the topological 
origin of the disclination defect. The elastic flux due to pentagonal apical disclination represented by abelian gauge 
field W is given by (see [11,15]) 

^W^dx^ = ^. (3) 

As is seen, the elastic flux for pentagonal defect is exactly 1/6. It was shown in [15] that the electronic eigenfunction 
acquires the additional Aharonov-Bohm-like phase due to nontrivial topology of elastic surface with a disclination 
vortex. In our case, this phase will appear both in the Bloch (1) and the Dirac (2) wavefunctions. 

To exclude the discontinuous boundary conditions for spinor fields on the sector sides (lines OB and OB' on Fig.l) 
let us introduce two additional fields. First, constructing a boundary condition for the spinor components on the line 
OB = OB' one has to take into account the exchange between A and B sublattices which takes place for the sector 
angle 7r/3. This can be done by using an appropriate boundary condition for the K spin part in the form (see [12]) 

iP{ip + 2w) = -TiP{ip), T = iT2 = exp{i^T2) (4) 

where the Pauli matrix T2 acts on the K part of the spinor components, and ip is the polar coordinate on the membrane 
{0 < ifi < 2tt). One can introduce the non- Abelian gauge field a to eliminate the exponential factor. In this case, the 
circulation of the field a is written as [12] 

/ 27r 

(b a^dx'^ = T2 — . (5) 

As is seen, this field adds ±1/4 to the total flux. It should be mentioned that the circulation of the gauge field a is 
governed by topology of the lattice and docs not depend upon geometry of the structure. Thus these operators and 
fluxes can be obtained solely from the lattice structure. 

The last field to be introduced is the frame rotation field Q which is equivalent to the spin connection (see [12]). 
This field (analog to the metrical connection coefficients) rotates the frame by the angle 27r/6, and its circulation 
along the contour A' A is determined by the following condition: 

jq^dx^ = -cj3^. (6) 

Notice that the spin connection does not contribute to the total flux. 



B. The covariant description and the Dirac equation on the curved surface 

So far we considered the problem on the plane by using the discontinuous planar coordinates (with the borders of 
the sector OB and OB'). Instead, one can also use the continuous coordinates on the Riemannian surface r(r, (/?), 
Q < (f <2ti (see Fig.l). To incorporate fcrmions on the curved backgroimd we need a set of orthonormal frames {ca} 
which yield the same metric, g^v, related to each other by the local 5*0(2) rotation, 

e^^e',=Kie0, A^ e 50(2). 

It then follows that = e^e^^Q^ where e'^ is the zweibein, with the orthonormal frame indices being a, /3 = {1, 2}, 
and coordinate indices fJb,y = {1,2}. As usual, to ensure that physical observables are independent of a particular 
choice of the zweibein flelds, a local so(2)-valued gauge fleld must be introduced. The gauge held of the local 
Lorcntz group is known as the spin connection. For the theory to be self-consistent, the zweibein flelds must be chosen 
to be covariantly constant [20], 
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which determines the spin connection coefficients explicitly 

[oj^r^ = e'^.D^e^r (7) 

Finally, the Dirac equation (2) on a surface S in presence of the U{1) external gauge field and the gauge field 
is written as 

n"e/(V^ - ia^ - iW^)i> = E^, (8) 

where = 9^ + il^ with 

i^/^ = ^'^"/[7a,7/3] (9) 

being the spin connection term in the spinor representation. Notice that = —iQ^ which justifies the above- 
mentioned relation between the frame rotation field and the spin connection. The spinor in (8) has the form = 
{F^ F4 " Fb )^ where F{r) are envelope functions, 

The matrix T2 acting in K-spin space appears in (8) only through a^. Therefore (8) can be easily diagonalized and 
we arrive at the two-component Dirac equations in the form 

-ia''e^{V^-ial~iW^)i^'' ^E^\ (10) 

As is seen, the coupled pair of equations (8) is reduced to the decoupled one describing "K-spin up" {K^) and "K-spin 
down" {K^) states. The field is determined by a condition 

j> aldx" = 

with the sign plus (minus) taken for k = (k = K^), respectively. Notice that after diagonalization the four- 
component spinor if) is found to be decomposed into "upper" and "lower" doublet components, and V^^> each 
of them transforms via SU(2). 

C. The continuum model for the icosahedral fullerene 

According to the Euler's theorem, the fullerene molecule consists of exactly twelve disclinations. Generally, it is 
difficult to take into account properly all the disclinations. There arc two ways to simplify the problem. First, one can 
consider a situation near a single defect (similar to [11]) taking into account that each defect in the fullerene can be 
simulated by two fluxes: K spin flux (5) and the elastic flux (3). In the case of sphere, however, the most appropriate 
approximation is to introduce the effective field replacing the fields of twelve disclinations by the field of the magnetic 
't Hooft-Polyakov monopole with a constant flux density and the half-integer charge A [10]. 

The total flux of the monopole Air A is equal to the sum of fluxes from all the disclinations. The procedure of 
summing up non-abelian fictitious fiuxes from apical defects placed at different points of the graphite cones was 
presented in [12] (so-called "n — m" combination rule). Namely, for the field a a linear integral circulating any even 
number of defects on the lattice is determined by {NTr/2 — (27r/3)[(n — m) (mod 3)]), where A'^ is a number of 
defects, n and m are numbers of steps in positive and negative directions, respectively (see [12]). The directions 
rotated by 27r/3 are considered to be identical. This gives a natural (n, m) classification of two-pentagon lattices: 
those for which n = m (mod 3), and those for which n ^ m (mod 3). This approach is suitable for the case of a 
sphere. 

As is shown below, the eigenf unctions of the low-energy levels oscillate not too fast with a distance, therefore the 
effective "monopole- like" approximation is valid for small quantum numbers (and near the Fermi energy). So, one 
can introduce the continuous field created by the 't Hooft-Polyakov monopole placed at the center of the sphere. In 
this case, the contour of integration is represented by two circles around the poles of the sphere as shown in Fig. 2. 
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FIG. 2. The integration path on the sphere. The path shown at the center encircles a few fivefolds. The "total" path of 
the integration, encircling all twelve pentagons, includes two contours around the poles: negatively directed for the north pole 
(shown) and positively directed for the south pole (not shown). 



Notice that for fullerenes with the full icosahedral symmetry (Ih) the combined flux turns out to be a sum of fluxes 
due to any pair of defects. Thus, the combined flux does not depend upon the arrangement of pentagons because 
any corresponding fragment of the lattice turns out to be of the above-mentioned class (1, 1). In other words, in 
icosahedral fullerenes (n — m) (mod 3) = due to the mirror symmetry of the lattice. Finally, the continuous fields 
take the form 

4 = 0,a^ = ±- cos6i. We = 0,W^ = -cosO. 

It should be noted that for the Dirac field in the external potential provided by a monopole, the effective "charge" A 
involves the "isospin" matrix r' (see [21]). This matrix can be also diagonalized in the Dirac equation thus giving the 
additional sign ± to the whole charge. What is important, the coordinate behavior of this field is the same as for the 
spin connection field (cf. [22,10]). This fact allows us to find the exact analytical solution of the problem. Another 
significant fact is the presence of "isospin" matrix in the momentum operator. Similarly to [10] 

= -z(V<^ - iAr^cosO) + y cos6 + AcosOt^ (11) 

and for half-integer A this operator takes integer eigenvalues j = 0, ±1, ±2, .... 

III. THE ELECTRONIC STATES OF THE FULLERENE 

In accordance with the results of the previous section, the total "charge" is written as A = ±(a^ + W^)/ cos6 = 
±1/2, ±5/2. Therefore, the Dirac operator in (8) takes the following form: 

- . / „ cot 9 , <^V , 

V = -la^ide + — -) - i^{d^ - lAcosO). 
I sm (7 

The substitution 

leads to the equations for Uj and Vj 

-i{de + [i - A] cot0 ± -^)v,{d) = EujiO), 
2 smu 

-i{de + cote- ^)nM = Evj{0). (12) 

The square of the Dirac operator reads 
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= -[a^ide + ^) + - Acose)f = 

z sm u 



1 



1 



I smf 



)t (9 



{c^zj + '^jA) + cot^ 9. 



sin 6 4 sin sin t 

Let us write the equation V'^'tjj = E^tp by using the appropriate substitution x = cosO. From (13) one obtains 



{j - Axf - ja,x + i + a^A 



1-x^ 



Uj{x) 

vj{x) 



Uj{x) 
Vj{x) 



(13) 



(14) 



Taking into account the asymptotic behavior of the spinor functions, one can use the substitution 

il-x)°'{l + xfujix 



{i-xy{i + xyvj{x) 



where 



1 


J-A-'- 


1 

' ^=2 


, 1 

J-A + - 




3+A+\ 


. 


3+A-\ 



Then the equation (14) for the function Uj takes the form 



1, 



(1 - x'^)dluj + (2(/3 - a) - 2(a + /3 + l)x)d^Uj + [-2a(3 -a- 13- ^{f - + ^- A) + - ^juj = 0. 

This is equivalent to the Jacobi equation 

(1 - x'')y" + {B-A-{A + B + 2)x)y' + A„y = 0, 

Xn=n{n + A + B+1), 



(15) 



(16) 



(17) 



with A = 2a, B = 2/3, and n being any non-negative integer. In view of (16) and (17) one gets the quantization 
condition 



A„ = n{n + 2{a + /?) + 1) = -2a(3 - a - f3 - ^{f - A^ + A) + E^ - ^. 



(18) 



Taking into account (15), one obtains the energy spectrum in the form E^ = {n + a + P + 1/2)^ — A^. When A = 0, 
one has a = (l/2)|j — 1/2|, /? = (l/2)|j + 1/2| and the energy spectrum is identical to the one found in [22] for the 
Riemann sphere without a monopole. 

In a similar manner one can study the equation for the function v. This gives another energy spectrum E^ = 
{n + -y + 5 + i)^ — A^. It should be mentioned that both solutions (for u and v) should satisfy (12). This is possible 
if the condition a + /3 = 7 + 5 holds true and, on the other hand, if one of the energies in the equations for u and v 
becomes zero. 

Let us consider the first case. The possible values of j are found to be \j\ > \\A\ + 1/2|. For E > one obtains the 
spectrum 



and the eigenfunctions 



El = {n+\j\ + l/2r-A^ 

uj = c^{i-xni+xfp^'^'^^, 
vj = c,ii-xni+xyp^^''^ 



(19) 



(20) 



P„ are the Jacobi polynomials. The unit of energy here is hVp/R where R is the fuUerene radius. One should note 
that (19) is not allowed for analysis of zero-mode states. In particular, the degeneracy of the zero-mode state can not 
be calculated using (19). According to (12), the factors C„ and in (20) are interrelated. Indeed, for j > one has 



- x^)d, - 27 - 2/3x]ap2^'2/3-i ^ ^^(^ ^ a;)C„p2^-i'2/3, 

- x^)d, + 2/3 - 27ar]C7„p27-i,2/3 = ^^(^ _ x)c,p2i,20-\ 



(21) 
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Setting a; = 1 in the first equation and using the definition 



one finally gets for j > 



J^:'''(l) = ^r(aVl|n? ' ^n = ±\/(27 + n)(2/3 + n), 



C„ = -^asigni^„ [ j + l/2-A ) [ n + j+A+1/2 ) 



In general case, for arbitrary signs of E and j, we obtain 



/I + 1/2- A sign j7 ^^+1^1+^+ 1/2 

In the second case, the possible values of j and n are determined by \j\ < \ \A\ — 1/2|, n = 0. One gets exactly one 
zero-mode at fixed j and positive fixed A 

uo = 0, vo= Cy{l - xy{l + x)\ (24) 

where the relation Pq'^'^^{x) = const is taken into account. Accordingly, if A < there exists only one zero-mode 
solution Uo- Thus, for all possible values of j and all possible positive values of A there exists exactly six different 
zero-mode solutions Vo- 
lt should be noted that both in (20) and in (24) the replacement j — j, A — » —A is equivalent to the exchange 
u — » z;, f — > M (cf. (12)). This means that "isospin up" and "isospin down" components of the spinor on a sphere 
turns out to be physically equivalent up to the redefinition of the quantum number j and a unitary transformation. 
Therefore, one can restrict consideration to either component, for instance, to "isospin up". 

From (19) and (24) one can calculate the energy spectrum. The possible "charges" arc A = —1/2,5/2, so that the 
first four levels arc (in units of hVp / R) the following: E = 0,1.41,2.45,3.46. Their degeneracies are g = 6,2,6,6, 
respectively. It is interesting to compare these results with tight-binding calculations. However, two preliminary 
remarks should be done. First, the continuum model is correct for the low-lying electronic states. Second, the validity 
of the effective field approximation for the description of big fuUerenes is not clear yet. In fact, the essence of this 
approximation is to take into account the isotropic part of long-range defect fields. Probably, for bigger fuUerenes 
both the anisotropic part of the long-range fields and the influence of the short-range fields due to single disclinations 
should be properly involved. Therefore, the exact values of the energy levels do not agree well with those presented 
in [4-7]. At the same time, we can surely verify both the existence of quasi-zero modes found for spherical fuUerenes 
in [4-7] and their 6-fold degeneracy. There is also a good qualitative agreement in observed scaling of the energy gap 
between the highest occupied and lowest unoccupied energy levels with the size of the fuUerene. 



IV. CONCLUSION 



In this paper, we have studied the electronic states of the icosahedral fuUcrcne within the continuum field-theory 
approach. The influence of the disclinations is taken into account by introducing an effective field due to magnetic 
monopole placed at the center of a sphere and having a total "charge" A- The flux due to monopole is a sum of two 
fluxes: (i) the K-spin flux and (ii) the elastic flux due to nontrivial topology of the surface with a disclination (the 
pricked out point on the surface). An exact analytical solution of the problem is found and the explicit form of the 
zero-energy modes as well as of the energy spectrum is presented. 

It should be noted that our approach differs from previous studies of fullerene molecules within the continuum 
models [9-11,23]. The effective monopole field introduced in [9] is identical with the K-spin field a^. Notice that a 
similar to [9] monopole fleld was used in the continuum model of the spheroidal fuUerenes [23]. Neither in [9,10] nor 
in [23] the analytical solution was found. The presented model differs also from [11] where the gauge field due to K 
spin was neglected and only one disclination on a sphere was described. 

Introducing the gauge field due to elastic vortex we obtain some principally new results. First of all, the effective 
monopole charge A takes two different values within the proposed model (—1/2 and 5/2) rather than ±3/2 for 
icosahedral fuUerenes in [9]. In turn, this finding affects the energy spectrum which is a combination of spectra for 
these two charges. Besides, the eigenfunctions are characterized by different from [10] conditions for the momentum 
3- 
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Notice that the more precise description of the fuUerenes requires inclusion of the electron-phonon interaction. In 
the simplest form, the role of this interaction was considered in [9]. It was shown that the energy levels become 
shifted and lose a symmetry around the Fermi level. The similar effect is expected in our model. An interesting 
open question is the electronic structure of other (non-Ih) types of spherical fullerenes as well as of non-spherical (e.g. 
elliptical) ones. Probably, the introduction of the monopole-like fields will be also a good approximation, at least for 
the low-energy states. 

We would like to acknowledge M. Pudlak and S. Sergeenkov for useful discussions and comments. This work has 
been supported by the Russian Foundation for Basic Research under grant No. 05-02-17721. 
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